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^ ABSTRACT 

The radio galaxy Fornax A (NGC 1316) is a prominent merger remnant in the outskirts 
of the Fornax cluster. Its giant radio lobes suggest the presence of a powerful AGN, 
and thus a central supermassive black hole (SMBH). Fornax A now seems to be in a 
transition state between active black hole growth and quiescence, as indicated by the 
strongly declined activity of the nucleus. Studying objects in this evolutionary phase 
is particularly important in order to understand the link between bulge formation and 
black hole growth, which is manifested in the M,-cr relation between black hole mass 
(— I and bulge velocity dispersion. So far a measurement of the SMBH mass has not been 

^ ^ possible in Fornax A, as it is enshrouded in dust which makes optical measurements 

impossible. We present high-resolution adaptive optics assisted integral-field data of 
Fornax A, taken with SINFONI at the Very Large Telescope in the K band, where the 
influence of dust is negligible. The achieved spatial resolution is 0.085 arcsec, which is 
^ about a fifth of the diameter of the expected sphere of influence of the black hole. The 

—I stellar kinematics was measured using the region around the CO bandheads at 2.3 |i.m. 

Fornax A does not rotate inside the inner ^ 3 arcsec. The velocity dispersion increases 
towards the centre. The weak AGN emission affects the stellar kinematics in the inner 
0.06 arcsec only. Beyond this radius, the stellar kinematics appears relaxed in the 
central regions. We use axisymmetric orbit models to determine the mass of the SMBH 
in the centre of Fornax A. The three-dimensional nature of our data provides the 
possibility to directly test the consistency of the data with axisymmetry by modelling 
each of the four quadrants separately. According to our dynamical models, consistent 
SMBH masses M, and dynamical Ks band mass-to-light ratios T are obtained for 
all quadrants, with (M.) = 1.3 x 10^ Mq (rms(M.) = 0.4 x 10^ M©) and (T> ^ 
0.68 (rms(T) — 0.03), confirming the assumption of axisymmetry. For the folded and 
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averaged data we find M, = l.^tol^ x 10^ M© and T = 0.65t!^ [5p {3a errors). Thus 



^ 1^ the bl ack-hole mass of Fornax A is consistent within the error with the lTremaine et aP 

Mt-a relation, but is a f actor ^ 4 smaller than expected from its bulge mass 



([2002 
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Marconi & Hunt ( 2003 1 relation 
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1 INTRODUCTION 

Studies of the dynamics of stars and gas in the centres of 
nearby galaxies with a massive bulge component have estab- 
lished the presence of supermassive black holes (SMBHs) in 
the lO'^-lO^ Mq range. The mass of the central SMBH is 



* Based on observations at the European Soutliern Observatory 
VLT (076.B-0457(A)) and arcliival ESO La Silla (66.C-0310(A)) 
and NASA/ESA Hubble Space Telescope data (GO Proposal 
7458), obtained from the ESO/ST-ECF Science Archive Facility 
t E-mail: nnowak@mpe.mpg.dc 



tightly correlated with the bulge mass or luminosity (e.g 
Marconi &: Hunt||2003 1 and with the bulge velocity disper- 



sion a (Gebhardt et al. 2000b Ferrarese & Merritt 20001 



These correlations suggest that bulge evolution and black 
hole growth are closely linked. Indeed there is increasing 
theoretical evidence that galaxy merging (or other processes 
that lead to gas inflow, like secular evolution), nuclear activ- 
ity and feedback are all somehow linked to bulge and SMBH 



evolution (e.g. 


Di Matteo et al.|2005 


Johansson et al.|2008 


Younger et al. 


2008 Hopkins et al. 


20081. To observation- 



ally constrain present theories of bulge and black hole evolu- 
tion, detailed studies of AGN and merger remnants in differ- 
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ent evolutionary stages are essential. NGC5128 (CenA) is 
presently the only galaxy with a powerful AGN that under- 
went a recent major merger which has a measured black-hole 



mass ( Silge et al.|2005 Marconi et al.|2d06 Neumayer et al, 
2007| 



A galaxy very similar to CenA is NGC 1316 (Fornax A). 
Fornax A is a giant elliptical galaxy located in the outskirts 
of the Fornax galaxy cluster. It is one of the brightest ra- 
dio galaxies in the sky with giant double radio lobes and 



5-shaped nuclear radio jets (Geldzahler & Fomalont||1984 1 



The peculiar morphology with numerous tidal tails, shells 
and loops ( |Schweizer||1980[ |1981[ ) suggests that a major 



merger happened about 3 Gyr ago ( Goudfrooij et al 



followed by some minor mergers (Mackie & Fabbiano 



20011, 



19981 



Its nucleus however is surprisingly faint in X-rays, which can 
only be explained if the nucleus became dormant during the 
last 0.1 Gyr (lyomoto et al. 19981. This makes it an ideal 



target for a dynamical black- hole mass measurement, as its 
stellar absorption lines are probably not or just marginally 
diluted by non-stellar emission from the AGN. 

Fornax A is classified as an intermediate form between 



core and power law galaxy in Lauer et al. ( 2007 1. The surface 
brightness approximately follows an r ^ -law (Caon et al 



|1994[ ). Large amounts of dust are present also in the in- 
ner few arcseconds, which affects optical spectroscopy. We 
therefore measure the two-dimensional stellar kinematics of 
Fornax A with the near-infrared integral-field spectrograph 
SINFONI at the Very Large Telescope (VLT) in the K band, 
where the dust obscuration is only about 7 per cent of the 
obscuration in the optical. We use this data to analyse the 



stellar kinematics in a way similar to Nowak et al. (20071 



Throughout this paper we adopt a distance to Fornax A 
of 18.6 Mpc based on HST measurements of Cepheid vari- 



ables in NGC 1365 (Madore et al. 19991. At this distance. 



1 arcsec corresponds to 90 pc. With a velocity dispersion 
of 226 km s~^ (mean a measured in an 8 arcsec aperture 
in §3.4) the estimated sphere of infiuence has a diame- 
ter of 0.46 arc sec and we would exp ect a SMBH mass of 
2.2 X 10** Mo ( [Tremaine et al.||2002 l. This is large enough 
to be resolved easily from the ground with adaptive optics. 

This paper is organized as follows: In §2 we present 
the data and the data reduction. The stellar kinematics of 
Fornax A is described in §3, the measurement of the near- 
infrared line strength indices in §4. The stellar dynamical 
modelling procedure and the results for the SMBH mass are 
presented in §5, and §6 summarises and discusses the results. 



2 DATA & DATA REDUCTION 
2.1 SINFONI data 

Fornax A was observed between 2005 October 10 and 12 as 



part of guaranteed time observations with SINFONI ( Eisen 



hauer et al. 2003 Bonnet et al. 20041, an adaptive-optics 



assisted integral-field spectrograph at the VLT UT4. The 
nucleus of Fornax A served as guide star for the AO correc- 
tion. We used the K band grating (1.95 — 2.45 ]xm) and, de- 
pending on the seeing conditions, the intermediate size field 
of view of 3 X 3 arcsec (0.05 x 0.1 arcsec spaxel"^, referred 
to as "lOOmas scale" in the following) or the high resolution 
mode with a 0.8 x 0.8 arcsec field of view (0.0125 x 0.025 arc- 
sec spaxel"^, "25mas scale"). The total on-source exposure 



time was 130 min in the highest resolution mode and 70 min 
using the intermediate plate scale, consisting of 10 min expo- 
sures taken in series of "object-sky-object" (O-S-0) cycles, 
dithered by up to 0.2 arcsec. In addition three O-S-0 cy- 
cles of 5 min exposure time each were taken using the low 
resolution mode (8x8 arcsec field of view, "250mas scale"). 
The SINFONI data reduction package SPRED 



( Schreiber et al. 2004 Abuter et al. 2006 1 was used to re- 



duce the data. It includes all common reduction steps nec- 
essary for near-infrared (near-IR) data plus routines to re- 
construct the three-dimensional datacubes. The data were 
skysubtracted, flatfielded, corrected for bad pixels and for 
distortion and then wavelength calibrated using a Ne/Ar 
lamp frame. The wavelength calibration was corrected us- 
ing night-sky lines if necessary. Then the datacubes were 
reconstructed and corrected for atmospheric absorption us- 
ing B stars that do not have strong spectral features in the 
region of the CO bandheads. As a final step all datacubes 
were averaged together to produce the final datacube. The 
data of the telluric stars were reduced likewise. The fiux cali- 
bration of the 250mas data was performed by comparison to 
the photometrically calibrated NTT/SOFI Ks band imag- 
ing (see §2.3 and Fig. |3|, that gives a Ks magnitude of 8.44 
in an 8 arcsec aperture. The lOOmas image was calibrated 
with respect to the 250mas image and the 25mas image with 
respect to the lOOmas image. Fig. [l] shows the resulting im- 
age of Fornax A in all three platescales, collapsed along the 
wavelength direction within the SOFI Ks band region. 

The spatial resolution, i.e. the FWHM of the point- 
spread function (PSF), is difficult to derive when dealing 
with high-resolution AO observations of diffuse objects. De- 
riving it from the science data itself can work quite well, if 
there are either stars in the field of view or the target is 
an AGN or starburst galaxy with certain spectral proper- 
ties (see e.g. Mueller Sanchez et al.|2d06 Davies et al.|2006 



|2007[ ). None of this applies to Fornax A. Therefore we de- 
rived the PSF by regularly taking an exposure of a nearby 
star with approximately the same R band magnitude and 
B — R colour as the central 3 arcsec of Fornax A. This pro- 
cedure is usually recommended to observers, but apart from 
being very time-consuming it can also be unreliable under 
certain conditions ( Davies et al.|2004 1 . The atmosphere may 
vary strongly with time, resulting in a measured PSF that 
is different from the real one. In addition the response of 
the wavefront sensor to the PSF star is, to some degree, dif- 
ferent from the response to the AO guide star (the nucleus 
in this case), because the nucleus is extended and therefore 
the fiux distribution on the wavefront sensor is different. 
We used this method as it was the only available option to 
measure the PSF. 

The high-resolution (25mas) data were taken under 
good and relatively stable ambient conditions with a near- 
IR seeing around 0.6 arcsec. The two-dimensional PSF (i.e. 
the image of the PSF reference star) is shown in Fig.|2^. The 
spatial resolution is FWHM « 0.085 arcsec. A Strehl ratio 
of ~ 45 per cent was reached. The reliability of this PSF is 
verified by a comparison of the 25mas SINFONI luminosity 
profile with an /fST/NICMOS F160W (camera NIC2, sam- 
pled at 0.075 arcsec px~^) luminosity profile. They agree 
well without further broadening of the NICMOS data (see 
below and Fig. |4]| . 

During the hours of degraded seeing (around ~ 0.9 arc- 
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sec on average in the near-IR) on October 12 we observed 
Fornax A using the intermediate plate scale (lOOmas). The 
spatial resolution inferred from the PSF is « 0.16 arcsec and 
the achieved Strehl ratio was ~ 18 per cent. The shape of the 
PSF (Fig. is rather asymmetric and may not represent 
the true PSF due to the changeable weather conditions. The 
seeing and the coherence time drastically improved during 
the observations with the lOOmas scale, so we switched to 
the 25mas scale for approximately one hour and only when 
the seeing and the coherence time then degraded again to 
similar values as before we could observe the lOOmas PSF. 
In order to get a better estimate of the real PSF we com- 
puted the kernel that transforms the 25mas image of For- 
nax A (binned 4 x 4 to match the lOOmas spaxel size) into 
the lOOmas image using the program of jGossl fc Riffeser] 
| 2002[ ). The 25mas PSF was then rebinned to the lOOmas 
spaxel size and convolved with this kernel. Kernel and re- 
sulting PSF are shown in Fig. [2j; and |2jl. The shape of the 
kernel and hence the convolved PSF are quite noisy. The 
reason is the low S/N of the 25mas image compared to the 
lOOmas image (cf. Fig. [T]). In order to find the closest match 
between the two images the kernel tries to incorporate the 
noise, too. Thus the resulting PSF has a S/N more similar 
to the 25mas image than to the image of a bright PSF star. 
Nevertheless the convolved PSF has about the same FWHM 
as the measured lOOmas PSF star. 

The PSF is an integral part of the dynamical modelling 
and it is therefore important to know its shape as accurate 
as possible. However, so far no studies have been performed 
that analyse systematically the dependence of the resulting 
M. and mass-to-light ratio on the PSF shape. Such a study 
is beyond the scope of this paper, nevertheless we will do 
the modelling, as far as lOOmas data are included, twice, 
using one time the measured PSF and the other time the 
convolved PSF. 



2.2 Longslit data 

Kinematics derived from longslit spectra are useful to con- 
strain the orbital structure at large radii, outside the small 
SINFONI field of view. Major-axis longslit data of Fornax A 
are available from different authors (ILonghctti ot al. 111998 



sure the higher-order Gauss-Hermitc coefficients /13 and 



(Gerhard 1993 van der Marel & Franx 1993) which quan- 



Arnaboldi et al.||1998| |Bedregal et 


al. 2006, Sagfia et al. 


2002 


I. Longhetti et al. ( 


1998 


1 and 


Arnaboldi et al. 


( 


19981 



tify asymmetric and symmetric deviations from a Gaussian 



velocity profile. Bedregal et al. ( 2006 1 measured all four pa- 



rameters from optical spectra, whereas Saglia et al. ( 2002 1 
used the CaT region. The velocities of all measurements are 
consistent, taking into account the different seeing values. 
The velocity dispersions in the central ~ 5 arcsec however 



that can be diverse. The authors used different correlation 
techniques, different spectral lines and slightly different po- 
sition angles. Fornax A contains dust, which might alter the 
kinematics measured from optical absorption lines. |Silge fc| 
Gebhardt ( 2003 1 found that optical dispersions tend to be 



larger than those measured from the CO bandheads, which 
could be an effect of dust. The largest dispersions are indeed 
those measured from optical spectra. We therefore used the 



kinematics of Saglia et al. (20021, as they measured all four 



parameters from the CaT line region. Their measurements 
agree best with our SINFONI measurements. They observed 
FornaxA in October 2001 at the Siding Spring 2.3m tele- 
scope using a 6.7 arcminx4 arcsec longslit and an expo- 
sure time of 60 min and determined the kinematics with the 
Fourier Correlation Quotient method ( [Bender et al.|[T994[ ). 
As they assume a position angle of 58° we adopt this value 
in the following analysis. The major axis profiles are shown 
in Fig. [17] and in [Beletsky et"aL| ( [2008| . 

2.3 Imaging 

To measure the SMBH mass in FornaxA, it is essential to 
determine the gravitational potential made up by the stellar 
component using photometric measurements with sufficient 
spatial resolution and radial extent. Therefore we combine 
high-resolution HST NICMOS F160W and ground-based 
wide-field Ks band imaging. The low-resolution Ks band im- 
age was taken with the near-IR imager/spectrometer SOFI 
on the 3.5m NTT telescope on La Silla, Chile. It has a 
field of view of 4.9 x 4.9 arcmin with 0.29 arcsec px~^ 
and a seeing of ~ 0.7 arcsec. It was dust-corrected with 
the method described in Appendix A using a J band SOFI 
image. The dust-corrected image is shown in Fig. [S] Note 
that the original SOFI image and the residuals are shown in 



measured only v and a from optical absorption lines (Mg6) 
and from the CaT region in the near-IR, they did not mea- 



Beletsky et al. (20081. We applied the same method to the 
HST NICMOS F160W image, which we corrected using a 
NICMOS FllOW image. The isophotal profiles of the dust- 
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a)25flniasPSF 




tii 25mBS PSF (rebinned lo lOOmas pixelsaale} 
convolved with kernel 




Figure 2. a) 25mas PSF derived by observing a star of about 
the same magnitude and colour as the AO guide star (i.e. the 
nucleus of Fornax A) during the observations, b) Same as (a) but 
for the lOOmas scale, c) Kernel that transforms the 25mas image 
of Fornax A to the lOOmas image, d) 25mas PSF from (a) rebinned 
to the lOOmas spaxel size and convolved with the kernel from (c). 
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Figure 3. The top panel shows the dust-corrected SOFI Ks band 
image of Fornax A (courtesy Y. Beletsky). The dust-corrected 
light distribution appears very smooth and regular and the 
isophotes are slightly boxy. For comparison the HST WFPC2 
F450W image (NASA/ESA HST data (GO Proposal 5990), ob- 
tained from the ESO/ST-ECF Science Archive Facility), which 
highlights the strong dust lanes, is presented in the bottom panel. 



corrected images were obtained following [Bender fc Moel^ 
[lenhofi] ( |1987[ ). The light distribution appears very smooth 
and regular and the isophotes are slightly boxy. The surface 
brightness profiles were then combined using the method of 
[Corsini et al.| ( |2008[ ) , matching the profiles in the region be- 
tween 1.5 arcsec and 5.0 arcsec and shifting the NICMOS 
profile by the computed amount (/xsofi — I^HST) '-'^ 
Ks band SOFI profile. The original and the combined pro- 
files are shown in Fig. |4] 

The combined photometry extends out to 80 arcsec and 
was deprojected for an inclination i = 90° using the program 
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width of the LOSVD bins for the given dataset (wavelength 
range, a, /13, /14, S/N, spectral resolution). Especially the se- 
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Figure 4. Surface brightness profiles of the HST NICMOS dust- 
corrected H band image (green crosses), the SOFI dust corrected 
Ks band image (red crosses) and the combination of the two sur- 
face brightness profiles (black line). The SINFONI 25mas K band 
image (scaled to match the HST NICMOS profile) is overplotted 
in blue. 



of Magorrian ( 1999 1 under the assumption that the galaxy 
is axisymmetric. No shape penalty was applied. The stellar 
mass density then can be modelled via p* = Tv, where v is 
the luminosity density obtained from the deprojection. The 
Ks band mass-to-light ratio T is assumed to be constant 
with radius. This assumption holds approximately in the 
central part of the galaxy for which we have kinematic data 
and where dark matter does not play a major role. 



3 STELLAR KINEMATICS 

The stellar kinematic information are extracted using the 



maximum penalised likelihood (MPL) technique of Geb- 
|hardt et al.| ( |2000a| ), which obtains non-parametric line-of- 
sight velocity distributions (LOSVDs) as follows: an initial 
binned velocity profile is convolved with a linear combina- 
tion of template spectra and the residuals of the resulting 
spectrum to the observed galaxy spectrum are calculated. 
The velocity profile is then changed successively and the 
weights of the templates are adjusted in order to optimise 
the fit to the observed spectrum by minimizing the func- 
tion Xp = X'^ + Q^T-', where a is the smoothing parameter 
that determines the level of regularisation, and the penalty 
function V is the integral of the square of the second deriva- 
tive of the LOSVD. We fitted only the first two bandheads 
^^CO(2 - 0) and ^^CO(3 - 1), i.e. the spectral range be- 
tween 2.275 \ira and 2.349 ii.m rest frame wavelength. The 
higher-order bandheads are strongly disturbed by residual 
atmospheric features. At wavelengths A < 2.29 \i.m most ab- 
sorption lines are weak such that extremely high S/N spec- 
tra would be needed. The strongest absorption line in that 
regime is NaT, but as its strength is much higher in early- 
type galaxies than in single stars due to enhanced silicon 
( Silva et al.||2008 l it cannot be fitted by the template spec- 
tra. 



lection of the smoothing parameter a is a crucial step ( Mer- 
ritt 1997 Joseph et al.| 20011. If a is chosen too high, the 
LOSVD is biased toward a flatter shape and if it is too small 
the LOSVD is too noisy. In order to find the best initial val- 
ues for the Fornax A data Monte Carlo simulations on a large 
set of model galaxy spectra have been performed. These are 
described in detail in Appendix B. 

Based on the simulations we conclude that with MPL 
reliable LOSVDs can be obtained from the first two CO 
bandheads of the SINFONI data of Fornax A when the cho- 
sen a is between 1 and 10 and the S/N > 30. The S/N of 
our data is very high (on average ~ 70 for the 25mas and 
the 250mas data, and ~ 140 for the lOOmas data). We are 
using a = 8 for the 25mas and 250mas data and a = 6 for 
the lOOmas data. 



3.2 Kinematic template stars 

We built up a small library of late-type (K and M) kinematic 
template stars which were observed with SINFONI in the K 
band during commissioning and GTO observations between 
2004 and 2006. In total we have nine kinematic templates for 
the 25mas scale (spectral resolving power R ~ 5000) , twelve 
templates for the lOOmas scale [R ~ 4500) and ten tem- 
plates for the 250mas scale {R ~ 4000). As shown in |Silge fc| 
Gebhardt ( 2003 1 a correct value of the velocity dispersion a 



can only be obtained from fits to the CO bandheads, if the 
template has about the same intrinsic CO equivalent width 
(EW) as the galaxy. The CO EWs of the Fornax A spectra 
were measured using the definition and the velocity disper- 
sion correction of |Silge fc Gebhardt] ( |2003[ ). The resulting 
values are between 13. OA and 14. 5A. 

After excluding all stellar templates with EWs far be- 
low or above the measured range, five templates remained 
for the 25mas scale, five for the lOOmas scale and four for 
the 250mas scale. As a cross-check and to avoid template 
mismatching the kinematics was extracted first with single 
templates. The results agreed, so no other template star had 
to be excluded from the sample. Tab.[T]shows the used kine- 
matic template stars with the according spectral types and 
CO EWs. For a better comparison with the line strength in- 
dices discussed in §4 the CO EWs were also measured using 
the definition of |Silva et aL] ( |2008[ ). 

3.3 Error estimation 

The uncertainties on the velocity profiles are estimated using 
Monte Carlo simulations ( Gebhardt et al.|2000a |. A galaxy 
spectrum is created by convolving the template spectrum 
with the measured LOSVD. Then 100 realisations of that 
initial galaxy spectrum are created by adding appropriate 
Gaussian noise. The LOSVDs of each realisation are deter- 
mined and used to specify the confidence intervals. 



3.1 Initial parameters 

Before being able to obtain the final LOSVDs it is neces- 
sary to quantify possible systematic offsets and to constrain 
the best initial values for the smoothing parameter and the 



3.4 The kinematics of Fornax A 

The SINFONI data are binned using a similar binning 



scheme as in Nowak et al. ( 2007 1 with five angular bins per 



quadrant and a number of radial bins (seven for the 25mas 
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Table 1. CO equivalents widths of the stellar kinematic tem- 
plate stars. 



Name 


spectral type 


CC 


d 


CC 


b 


scale 


HD12642 


K5 III 


12 


7 


17.09 




25 


HD163755 


K5/M0 III 


15 


6 


21.27 


25 


HD179323 


K2 III 


11 


1 


14.08 




25 


HD198357 


K3 III 


10 


8 


13.^ 






25 


HD75022 


K2/3 III 


12 





15.06 


25, 


100 


HD141665 


M5 III 


15 


8 


22.26 


100, 


250 


HD181109 


K5/M0 III 


12 


9 


17.79 


100, 


250 


HD201901 


K3 III 


11 


4 


14.43 


100, 


250 


SA112-059^ 


MO III 


12 


8 


17.61 


100, 


250 



Note. The CO equivalent widths are given in units of A. The 
spectral types are taken from Wright ct al.| | |2003[ l where avail- 
able. In the last column the SINFONI platescales, in which the 
stars have been observed, are given. 



" Definition of'Silg( 
b 



Definition of ^Silva et al ^ 
N. Neumayer, priv. comm 



Gebhardtl |2003 I 
2008(1 






(b) . 


ti„,i+\/ 





2.30 2.32 
X(nm) 



-500 500 
V (km/s) 



Figure 5. a) Example fit to a 25mas spectrum and b) the ac- 
cording LOSVD with error bars. 



scale, ten for the lOOmas scale and 13 for the 250mas scale). 
The centres of the angular bins are at latitudes i9 = 5.8°, 
17.6°, 30.2°, 45.0° and 71.6°. The stellar kinematics then 
was extracted for all scales using the appropriate templates 
and Q. Fig. [5] shows as an example the fit to a 25mas spec- 
trum and the according LOSVD. 

For the dynamical modelling we use these LOSVDs and 
no parametrized moments, but for illustration purposes we 
fitted the LOSVDs with Gauss-Hermite polynomials and we 
show the 2D fields of t;, cr and the Gauss-Hermite coefficients 
/13 and hi in Fig.[6]for all platescales. The amplitudes of the 
errors are on average about 8 km s~^ {v and a) or 0.025 
(/13 and hi) for the 25mas data, 4 km s~^ or 0.015 for the 
lOOmas data and 9.5 km s~^ or 0.033 for the 250mas data. 

On large scales (250mas, Fig. ^) the galaxy is clearly 
rotating and the mean velocity dispersion is high (a ~ 
226 km s~^), while in the central ~ 3 arcsec (lOOmas, Fig. 
^j^) no clear rotation is visible and the dispersion is lower 
(cr « 221 km s"^). The 25mas velocity field (Fig. [eji) is ir- 
regular and with a peak in the central bins. The velocity 
dispersion also rises towards the centre, but the maximum 
does not coincide with the photometric centre. It is located 
~ 0.05 — 0.1 arcsec south of the centre. North of the centre 
there is a second but less strong maximum. The two maxima 
are separated by a narrow, elongated cr-minimum. Outside 



the central region the mean dispersion of the 25mas field is 
lower than in the larger fields (cr ~ 218 km s~^). /13 and /14 
are on average small and positive. 

In addition to the difficulty in explaining the structures 
of V and a in the centre, it is also surprising that neither 
the V nor the a peak can be seen in the lOOmas data. We 
investigated if this is due to resolution in the following way: 
The 25mas data were binned 4 X 4 to get the same spaxel 
sizes as in the lOOmas scale. Then the binned data were con- 
volved with the kernel (Fig. ^) that transforms the 25mas 
SINFONI image of Fornax A into the lOOmas image, before 
being binned according to the radial and angular scheme 
described above. Fig. [7|i shows the resulting 2D kinematic 
maps compared to the inner 0.8 arcsec of the lOOmas data 
(Fig. [7]d). The structures seen in the 25mas data mostly 
disappeared after the convolution and the convolved 25mas 
kinematics now look very similar to the lOOmas kinematics. 

The interpretation of the features seen in the uncon- 
volved 25mas data is more difficult. In the following possible 
explanations are discussed. 

3.4-1 Technical aspects 

The features detected in the 25mas data, especially the cen- 
tral velocity rise, are not due to template mismatching, as 
they remain also when other template combinations or sin- 
gle templates are used. Of the nine 25mas templates, two (a 
KIV and a K3V star with CO EWs < 5A) cannot fit the 
data at all. With each of the other templates the kinematics 
of Fig. [6|i can be reproduced at least qualitatively indepen- 
dently of the specific choice of the smoothing parameter. 

The results can also be reproduced when fitting slightly 
different wavelength regions, as long as the CO bandheads 
are included. Measuring the kinematics from the other ab- 
sorption lines (Cal, Fel and Mgl) does not work because 
of the weakness of the lines compared to CO. Here a much 
higher S/N would be required. 

The central bins are very small and only cover a limited 
number of spaxels. Therefore one or two spaxels alone could 
cause the increase of v, because e.g. of an incorrect sky sub- 
traction or bad pixel correction in the CO bandhead region 
of these spaxels. In addition, as the S/N of the single spax- 
els decreases relatively slowly from the centre to the outer 
parts, the S/N of the central bins is lower (~ 50) than at 
larger radii (60 — 80), so the errors are larger and the large 
velocity could be just a statistical fluctuation which would 
occur at larger radii as well if the S/N was equally low. A 
single-spaxel effect can be in principle verified or excluded 
by simply fitting the unbinned spectra. The problem is, that 
the S/N is around or below 30 which does not allow a re- 
liable measurement of the kinematics. We first used MPL 
with a large a and then, as a consistency check, fitted sim- 
ple Gaussians to the unbinned spectra. A central v peak is 
present in both cases. 

3.4-2 A central disc or star cluster 

Central velocity dispersion drops are common in spiral 
galaxies (e.g. Peletier et al. [2007) and usually associated 



with central discs. Central discs can be formed from gas in- 



fiow towards the centre and subsequent star formation ( Woz- 



niak et al. 20031. In early-type galaxies cr-drops are rarely 
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Figure 6. a) v, a, hs and hi of Fornax A in the central 0.8 arcsec (25mas scale). The major and the minor axis as well as the 
numeration of the different quadrants are indicated in the velocity field (left), b) shows the same as (a) but for the lOOmas scale. 
The 25mas field of view is marked for a better orientation, c) shows the same as (a) but for the 250mas scale. Both the 25mas 
and the lOOmas field of view are marked. 
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(a) 25mas binned to lOOmas spaxel size and convolved with kernel: 




(b) lOOmas, inner 0.8": 
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Figure 7. The 25mas data convolved with the kernel from Fig.[2|; (upper row, (a)), compared to the central part of the lOOmas 
kinematics (lower row, (b)). 



found. NGC 1399 is the only case where the observed central 
a-drop is discussed in detail fHoughton et al. 2006 Gebhardt 
|et al |2007| [Lyubenova etal. 2008). Lyubenova et al.| ( |2008| ) 
conclude that it is a dynamically distinct subsystem with a 
different stellar population. 

Fornax A is not a clearly elliptical galaxy, but can be 
classified as a "peculiar" galaxy, as some remaining features 
like shells and ripples from a major gas-rich merger about 
3 Gyr ago are still clearly visible (Schweizer 1980 



1981 



[Goudfrooij eraL]p00T| |. There are also hints that a minor 
merger with another gas-rich galaxy happened ~ 0.5 Gyr 
ago ( Mackie & Fabbiano 1998") . Thus it would not be surpris- 



ing if remainders of the merging processes were also found 
in the centre, like a disc produced from merger-triggered gas 
inflow or an infalling disc- or cluster-like central part of one 
of the merging components. A nuclear disc could explain the 
relatively thin (< 8 pc) and apparently elongated structure 
of the a-drop. With a width of only < 8 pc this is probably 
the thinnest cr-drop ever detected. It is even narrower than 
the drop in NGC 1399 and also far less extended than the 
typical drops in spiral galaxies. As its size is comparable to 
the spatial resolution it is possibly unresolved. On the other 
hand rotation would be seen if there really was a disc. 

The centre could still be unrelaxed and in the process of 
merging. This could induce signatures in the kinematics like 
those seen in Fornax A. Simulations of unequal- mass merg- 
ers with a SMBH in the primary galaxy and a significantly 
smaller secondary galaxy (without black hole) were per- 
formed by Holley-Bockelmann & Richstone ( 2000 1 . Based 



on this it depends critically on the orbital decay trajectory 
of the secondary, and also somewhat on the black-hole mass 
and on the mass ratio of the galaxies, whether the secondary 
galaxy is destroyed by the merger or not. The time-scale for 
a stellar cluster falling in on a purely radial orbit - also 
suggested by Gebhardt et al. (20071 as one possible expla- 
nation for the (T-drop in NGC 1399 - however seems to be 
much shorter (of the order of about 10^ yr) than the date of 
the last merger (~ 10^ yr), if the conditions in Fornax A are 
similar as in Holley-Bockelmann & Richstone (20001. For 



nonradial mergers the formation of a rotating stellar disc 
is possible. As we do not know any of the details of the 
merger history with desirable precision for Fornax A, and as 
the time-scales seem to vary depending on the initial pa- 
rameters of the merger simulation, it is impossible to tell 
whether the centre is still unrelaxed or not. 



3.4.3 The effect of dust 

Fornax A contains a lot of dust features (see lower panel of 
Fig. [3| and this dust might alter the kinematics. On the 
other hand observations in the K band should minimize the 
effects of dust and indeed already the uncorrected SOFT Ks 
band image shows only very little amounts of dust compared 
to optical HST images. Some dust is located close to the cen- 
but not in the central 



tre, but not in the central ~ 2 arcsec. Shaya et al. ( 1996 \ 
find a low colour excess in the central pixels of an HST 
WFPC V — I image and conclude that here the dust extinc- 
tion is only very light. Consistently the SINFONl 25mas and 
lOOmas images also do not show any signs of dust (cf. Fig. 
[l|. The 250mas image shows hints of dust in the SE region, 
spatially coincident with the high-a region just outside the 
lOOmas field of view. Therefore the dust outside the central 
2 — 3 arcsec could be the reason for the larger average a of 
the 250mas data, in agreement with the findings of [Silge fc] 



Gebhardt (20031 



Due to the lack of dust in the 25mas image dust does 
not seem to be a plausible explanation for the central a 
structure. However the scale of the effect of dust on the 
kinematics might depend on the location of the dust along 
the line of sight. E.g. a certain amount of dust directly be- 
hind the centre could have a noticeable effect without being 
noticed in the photometry at the given S/N. To further sub- 
stantiate this we tested the effects on the kinematics when 
dust obscures light from stars behind or in front of the centre 
of a galaxy using an N-body model realisation of the best- 
fitting dynamical model (cf. §5.2.3) as described in Thomas 



et al. (20071. We calculated the kinematics of the model by 



© 2008 RAS, MNRAS 000,[T|f23l 



The SMBH of Fornax A 9 



obscuring the central part at different radii and with differ- 
ent opacities and found that only a noticeable decrease in 
0" can be induced by putting a very thick layer of dust in 
front of the centre, but in this case the surface brightness 
would be reduced by more than a magnitude. Therefore we 
can exclude that the structure of the 25mas kinematic fields 
is produced by dust. 



3.4-4 Stellar populations 

Another explanation might be that the stellar population 
is different in the central region. Blue stars in the nucleus 
could in principle also be responsible for the low colour ex- 

The K band absorption line 
4 in de- 



Shaya et al. 



(1996| find 

equivalent widths for Fornax A are presented in 
tail. The lOOmas maps of the absorption line indices (Fig. 
[8| clearly show that the line indices change within the cen- 
tral ~ 0.4 arcsec. All indices show a decrease in this region, 
CO as the deepest absorption feature shows it most clearly. 
Unfortunately this trend is not as clear in the 25mas indices, 
as here the S/N is only about half as in the lOOmas data and 
the scatter therefore is much larger. 

Thus a dynamically colder subsystem with accordingly 
different line indices could be an explanation of the distorted 
kinematics, but the S/N of the 25mas data unfortunately is 
not large enough to resolve any spatial coincidence with the 
V peak or a drop. Similar conclusions have been made by 
Lyubenova et al.| ( |2008[ ) for NGC 1399, where the central 
cr-drop coincides with a drop of the NaT and CO indices. 

If a changing stellar population in the central region has 
an effect on the kinematics or not was tested with the N- 
body model mentioned in the previous section. Making the 
stars in the central 0.1 arcsec ten times brighter (correspond- 
ing to a ten times lower mass-to-light ratio) results only in 
very small changes of Act — 4 km s~^ and A/14 = —0.002, 
both within the errors of the data and much smaller than the 
observed variations. 100 times brighter stars would produce 
noticeable changes of Act — 25 km s~^ and A/i4 — —0.015. 
Making the stars darker does not have an effect at all on the 
kinematics. As there is no photometric evidence of a very 
bright star cluster in the centre, we conclude that a stellar 
population alone, which is different from the surrounding 
population, is not the reason for the disturbed central kine- 
matics. 



3.4-5 AGN emission 

A decrease of the line indices would also be expected for 
AGN, where the absorption features are diluted by the non- 
stellar emission of the active nucleus. In this case also the 
measurement of v and a could be affected. As shown in 



Davies et al. ( 2007 1 the PSF can be reconstructed from the 
decrease of the CO equivalent width at the position of the 
point source. Comparing the PSF measured this way with 
the measured or the noisy reconstructed PSFs from Fig. 
[2] is difficult, because (1) we would need to measure the 
CO index of the unbinned spectra, which have a lower S/N 
than the binned ones and thus errors larger than the small 
~ 1.0 A decrease observed in the binned spectra, and (2) for 
the higher S/N lOOmas exposure we do not exactly know 
the true PSF. From the binned spectra we cannot deduce 
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Figure 9. H2 emission of Fornax A (lOOmas scale). The photo- 
metric centre is marked by a cross. 



the two-dimensional PSF. A rough estimate of the FWHM 
of the CO dip can be obtained from fitting the averaged CO 
indices of a ring (red points in Fig.[8|. The resulting FWHM 
of 0.17 arcsec in the lOOmas scale is, considering the large 
errors introduced by the binning, well in agreement with the 
measured PSF. This would favour the AGN scenario. 



It was found by lyomoto et al. (19981 that the AGN 



in Fornax A, which powered the very luminous radio lobes, 
became dormant during the last 0.1 Gyr. Still the nucleus 
shows weak radio (Geldzahler &: Fomalont||1984 1 and weak 
X-ray emission with a spectrum typical for a low-luminosity 
AGN (|Kim fc Fabbiano|2003[ ) and is bright in the UV ( [FaF 



biano et al.|1994" l. In high-S/N optical spectra several emis 



sion lines are present ([OIII], H/3, NI), but very weak (Beu- 
ing et al.||2002). In the SINFONI spectra there are no ap- 



parent emission lines typical for an AGN present at the 
given S/N. Only weak molecular hydrogen (1 — 0S(1) H2 
at A = 2.122 (xm) with a S/N around 5 is found in a re- 
gion north-east of the nucleus (see Fig. [9]l, but does not 
seem to be associated with the central source. The total H2 
flux in the lOOmas field of view is 5.65 x 10~^** W m~^ and 
in a 3 arcsec aperture centred on the continuum peak it is 
3.64 X 10"^* W m"^. The H2 could be excited by X-rays, as 
a similar northeast-southwest elongation is present in Chan- 
dra data. This X-ray emission was associated with hot inter- 
stellar medium with a temperature of kT = 0.62 keV by |Kim| 
|fc Fabbiano] ( [20031 ). 

A weak AGN present in Fornax A would also become 
noticeable by a slight change of the continuum slope, making 
the nucleus redder. This is due to the AGN UV-continuum 
which heats surrounding dust that emits in the near-IR 
(Oliva et al. 1999'). We indeed find that the slope of the 
continuum changes slightly in the inner ~ 0.1 arcsec by no 
more than 2 — 3 per cent. 

Some AGN show a certain coronal emission line in the 
region of the CO bandheads, which alters the measured kine- 
matics: [Ca VIII] at 2.3213 nm ( [PortiUa et al.|20 08: Davies] 
et al.|[2006| . As it is located directly at the left edge of the 
second CO bandhead, already a very small contribution can 
have a significant effect on the measured kinematics. We 
tested this using two different approaches. First we tried 
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Figure 8. Nal, Cal, <FeI> and CO indices of FornaxA in the central 1.5 arcsec (lOOmas scale). Tlie upper row sliows the 2D 
index fields, below the position-i diagrams are shown (where i is one of the four indices). The individual values are plotted in 
grey, while overplotted in red are the mean indices and their rms of the bins belonging to a ring of radius r. 



to measure the kinematics using only the first bandhead. 
The velocity map in this case has a much larger scatter and 
shows many peaks of 20-30 km s~^ in the field of view, a 
remains largely unchanged. When the first together with 
the third and fourth bandhead are used, the velocity map is 
smooth and shows no increase in the centre. The a map is 
qualitatively similar to the one measured with the first two 
bandheads, but the absolute value of a is much larger. The 
third and fourth bandheads are significantly distorted by sky 
emission lines and therefore do not give reliable results. The 
second approach was to use a standard star, convolve it with 
aLOSVD {v = Okms"\ a = 230 km s~\ /13 = ^4 = 0) and 
add the [Ca VIII] line with a typical FWHM of 150 km s^^ 
The results are shown in Figs. |10| and |ll[ A [Ca VIII] contri- 
bution of < 4 per cent would not be seen in the spectra or the 
fit at the given S /N, but can increase v hy < 25 — 30 km s~^ 
and decrease f by < 10 km s~^ {a measured with zero 
[CaVIII] contribution is already a few km s~^ lower than 
the real a, see e.g. Fig. |Bl| in Appendix B). Therefore a weak 
AGN with some [CaVIII] contribution is the most logical 
explanation for the observed velocity increase. The observed 
a decrease probably has for the most part a different origin 
(most probably a stellar population effect), as it is signifi- 
cantly larger than the < 10 km s~^ estimated from the sim- 
ulations. In addition the high-u region only partly coincides 
with the low-cr region. The hint of a [Ca VIII] line signature 
appears when stacking all spectra within the region covered 
by the innermost bins and comparing it to a combination of 
all spectra within a ring further out (0.11 < r < 0.17 arcsec; 
see Fig. [12 I. 



As accurate kinematics for the centre cannot be ob- 
tained from only one bandhead or by including the third 
and/or fourth bandhead, we will exclude the central bins in 
some of the the dynamical models. As the black hole sphere 
of influence is large, this solution will not degrade the relia- 
bility of the resulting black-hole mass significantly. 



4 LINE INDICES 

The optical line strengths of FornaxA have been analysed 
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Figure 12. Combined spectra of the central region (black) at 
radii r < 0.03 arcsec and of a ring at larger radii 0.11 < r < 
0.17 arcsec (green). In red the residuals between the two spectra 
are plotted, shifted by +1. The position where the [CaVIII] line 
is assumed to be is marked with an arrow. 



optical index gradients slightly decrease with radius. Un- 
fortunately their data have a resolution of only 1 arcsec, 
therefore they do not resolve the structure in the central 
0.4 arcsec. The age of a number of globular clusters in For- 
naxA was found to be about 3 Gyr ( Goudfrooij et al.|2001 1. 
Thus a plausible scenario is that the globular clusters formed 
during the last major merger event 3 Gyr ago (Goudfrooij 



et al. 2001 1 and that the young stellar population formed 



from infalling molecular gas also as a result of this merger 
(Horellou et al. 2001). 



The spectral features in the K band can in principle be 
likewise interpreted in terms of age and metallicity. Unfortu- 
nately there is not such a sophisticated theoretical spectral 
synthesis model available yet as for optical indices and stud- 
ies of K band indices are rare and often include only small 



samples. A recent study of Silva et al. ( 2008 1 investigates the 



extensively by Kuntschner ( 1998 2000 1. They found a young 



(around 2 Gyr) and very metal-rich stellar population. The 



behaviour of the K band indices of cluster stars of known 
age and metallicity and a sample of Fornax galaxies. Their 
published line indices of Fornax A are derived by extracting 
the mean spectrum over an aperture of 2 x 3 arcsec from our 
SINFONI lOOmas data. 

In Fig. [S] we show the two-dimensional K band line 
index maps (Nal, Cal, <FeI> and CO) of FornaxA for 
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Figure 10. Mean v, a, /13 and hn as a function of [CaVIII] contribution. Tiie error bars are the 67 per cent intervals obtained 
from fitting 100 realisations of the broadened template spectrum with added noise such that the S/N is similar to the central 
Fornax A spectra. 
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Figure 11. Fit (red) to broadened template star spectra (black) with different [CaVIII] contributions compared to the fit to the 
central 25mas Fornax A spectrum with the largest measured velocity (leftmost plot). 



Table 2. Mean near-IR line indices of Fornax A. 



this work 


Silva et al. 


1 2008 1 


Nal 

Cal 

FelA 

FelB 

CO 


4.61 (0.15) 
2.14 (0.21) 
1.44 (0.11) 
0.73 (0.09) 
17.61 (0.33) 


4.70 (0.14) 
2.48 (0.14) 
1.491 (0.066) 
0.843 (0.065) 
17.53 (0.30) 



Note. The indices are given in units of A. The middle column 
gives the near-IR line indices determined by averaging the 
values shown in Fig. [s] (upper row) and the according rms 
errors. The indices determined by [Silva et al.| l [2008[ l are given in 
the right column with their errors. 



indices show the same behaviour: a more or less clear dip in 
the centre and further out a negative gradient. This is most 
obvious in CO and NaT. Negative gradients are also present 
in the optical indices and it is therefore reasonable to assume 
that the K band index gradients can be likewise interpreted 
as age or metallicity gradients. |Kuntschner| ( |1998[ ) find for 
the Fornax ellipticals that the index gradients are linked 
to a metallicity gradient, and not an age gradient. They 
assume that the negative gradients in Fornax A are likewise 
due to a metallicity gradient only, but their data do not 
go very far out in radius and therefore they cannot draw 
secure conclusions. [Marmol-Queralto et al.| ( |2008[ ) quantify 
a metallicity dependence of the CO index, which supports 
this interpretation. 



the lOOmas scale, because the S/N is largest for this scale. 
The index definitions were taken from table 6 in lSilva et al.l 
( [2OO8I ). Instead of Fe I A and Fe I B we use <Fe I>=(Fe I A-|- 
FeIB)/2, as they are highly correlated. The definition of the 
CO index is different from the definition of the CO equiv- 
alent width used above ( jSilge fc Gebhardt||2003[ ). Here the 
continuum is defined by four blue continuum bands instead 
of a blue and a red one. The continuum on the red side is 
diminished by the CO absorption bands, therefore an ac- 
curate measurement of the CO index is possible by using 
blue continua only. Using the definition of |Silge fc Gebhardt| 
( |2003[ ) nevertheless is sufficient to decide quickly which tem- 
plates to use and to notice any particularities (see Table[l]for 
a comparison between the CO EWs of the kinematic tem- 
plate stars derived by the two different index definitions). 



As in |Silva et al.| ( |2008[ ) we broadened the SINFONI spectra 
to match the ISAAC resolution of their data and to be able 
to directly compare our results to their findings. 

The measured indices, averaged over the entire field of 
view, are in good agreement with the measurements of |Silva| 
et al. (20081, as shown in Tab. [2] In two dimensions, all 



What causes the dip in the very centre is less obvious. 
It might be one of the scenarios discussed above - AGN 
emission, dust or a different stellar population like a nu- 
clear cluster or a population of younger, bluer stars of a 



post-starburst population. Fabbiano et al. ( 1994 1 detected 



a UV-bright unresolved central source with r < 3 pc, which 
according to them could be explained e.g. by a cluster of O 
stars or an AGN. The FWHM of the CO dip is very similar 
to the spatial resolution and thus unresolved, corresponding 
to an upper limit in size of r < 4 pc. Thus it is likely that the 
UV emission and the CO dip are caused by the same struc- 
ture. Both early-type stars - due to the lack of typical late- 
type star absorption lines - and non-thermal emission would 
cause the observed decrease of all indices. Br 7 absorption 
is the only strong feature in early-type K band stellar spec- 
tra, Br 7 emission is observed in many galaxies with nuclear 
activity or star formation, thus this feature could be used 
to distinguish between these two possibilities. In the central 
SINFONI spectra we do not detect Br 7 neither in emission 
nor in absorption. Thus emission or absorption is either not 
present, too weak to be detected, or Br 7 emission and ab- 
sorption just compensate. [CaVIII] emission caused by the 
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AGN is also very well hidden (see Fig. 12 and discussion in 



§3.4.5). Due to the broader PSF its influence is eliminated 
in the lOOmas data. Other emission lines typically present in 
AGN spectra are absent as well in the central region. Weak 
H2 emission is present in the north-eastern region, but not 
correlated to the central source (see Fig. |9|. In order to re- 
solve the reasons that cause the central line index dip and 
kinematics, very high S/N observations with a better spatial 
resolution as well as self-consistent spectral synthesis models 
for the interpretation of the K band indices are needed. 



5 DYNAMICAL MODELS 



Based on the M,-a relation of |Tremaine et aL] ( |2002[ ) and our 

a measurements a black hole with a mass around 2 x 10* Mq 



would be expected. Earlier models by Shaya et al. ( 1996 1 



suggest a ten times higher black-hole mass, but they rely on 
larger a values. Davies (200011 observed Fornax A with the 



MPE 3D IFS (Weitzel et al 



1996 1 and measured v and a 



using the CO bandheads at 2.2 \im. Their a is lower than 



the values used by Shaya et al. ( 1996 1 and comparable to 



our values. Together with the models of Shaya et al. ( 1996 1 
they obtain a SMBH mass of < 10^ M©. However, they 
cannot exclude that radial anisotropy may account for this 
large mass. A lower limit of the central mass concentration 
of 75 M0 was estimated by |Fabbiano et aL] ( |1994[ ) based on 
the measured nuclear UV emission. 

We used the axisymmetr ic code of [Gebhardt et al.| 
( 2000a| pOOSl in the version oflThomas et al 1 1|2004|) to de- 



termine the mass of the SMBH in Fornax A. It is based on 



the Schwarzschild ( 1979 1 orbit superposition technique and 
comprises following steps: (1) Calculation of the gravita- 
tional potential of the galaxy from the stellar mass density 
p, using a trial black-hole mass M, and mass-to-light ra- 
tio T. (2) Generation of an orbit library for this potential 
and construction of a weighted superposition of orbits that 
best matches the observational constraints. (3) Repetition 
of the first two steps with different values for T and Al, 
until the eligible parameter space is systematically sampled. 
The best-fitting parameters then follow from a x^-analysis. 
Our orbit libraries contain around 2 x 7000 orbits. The de- 
projected luminosity density is a boundary condition and 
hence exactly reproduced, while the LOSVDs are fitted in 
25 velocity bins between —880 ... -I- 880 km s~^ with a bin 
width of ~ 75 km s~^. Special care was taken when imple- 
menting the PSF. Due to its asymmetric shape the PSF was 
not fitted, but the models were rather convolved with the 
two-dimensional image of the PSF reference star directly (cf. 
Nowak et al. 20071. All modelling was done with minimal 



regularisation. 

A big advantage of integral-field data compared to 
longslit data only is the assurance of whether the assumption 
of axisymmetry is legitimate by comparing both the kine- 
matics and the results of the dynamical modelling of all four 
quadrants. In case all quadrants produce the same results, 
the data can be folded and averaged, so that the errors are 
reduced. This results in a very large number of models. For 
Fornax A we calculated around 5000 models in total, cover- 
ing a broad parameter space in M, (0 ... 5 x 10* Mq) and 
T (0.55 . . . 0.95) for each quadrant, several combinations of 
datasets and PSFs using an inclination of 90°. At the begin- 




0.7 0.8 
T (Ksband) 

Figure 13. a-s a function of mass-to-light ratio T for models 
calculated with the longslit data alone, truncated at different radii 
J'trunc and without black hole. The models with the smallest 
are marked with a star. 



ning we computed some test models using either the longslit 
data only (§5.1.1) or the SINFONI data only (§5.1.2) in or- 
der to constrain the T range. As a second step each of the 
four quadrants was modelled separately using the longslit 
data and either the 25mas data (§5.2.1) or the lOOmas data 
(§5.2.2) or a combination of the two (§5.2.3). This step is 
essential as it reveals possible inconsistencies between the 
datasets, non-axisymmetries and other problems. The as- 
sumed diameter of the sphere of infiuence is around 0.46 arc- 
sec, resolved by both SINFONI datasets. The 250mas data 
were not included in the modelling, as they do neither re- 
solve the sphere of influence nor do they add signiflcantly 
more information in the outer parts that would justify the 
enhanced amount of computing time. After a careful exam- 
ination of the results and verifying the consistency of the 
data with axisymmetry, the four quadrants were averaged 
and modelled for the same combination of datasets as for 
the single quadrants. 



5.1 The stellar dynamical Ks band mass-to-light 
ratio T 

In order to check if the stellar mass-to-light ratios T ob- 
tained for individual datasets, agree and if they are consis- 



tent with T from population synthesis models (Maraston 
|1998[ p005| l, we determine T by modelling single datasets. 

5.1.1 T from longslit data 

First a number of models were calculated using exclusively 
the longslit data without black hole covering a broad range 
in T (in units of Mq /Lq^). The longslit data were truncated 
at different radii rtrunc between 14 arcsec and 45 arcsec. The 



resulting best-fitting T (Fig. 13 1 is roughly constant with 
rtrunc and « 0.75—0.8 for longslit data at rtrunc < 31 arcsec« 
2.8 kpc, which roughly correspon ds to the effective radi us of 
FornaxA {Rc = S6tli,2 arcsec, [Bedregal et alT||2006| . At 
larger truncation radii T increases strongly, a sign that here 
the dark halo starts to have an effect. This is consistent 



with the results of Thomas et al. (20071, who found dark 
matter fractions of ~ 30 per cent at 3 kpc in similar bright 
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Figure 14. Axq = — Xmin ^ function of T {Ks band) for 
models calculated with the SINFONI lOOmas data alone and a 
small number of black- hole masses M, . The coloured regions are 
the 1 . . . 5(T errors in Ax§. Each black dot represents a model. 



Coma galaxies. Thus in the following only longslit data at 
r < 31 arcsec will be used for the dynamical modelling. 



5.1.2 T from SINFONI data 

For the SINFONI lOOmas data alone the best-fitting T is 

Ks 



141, 



slightly lower (« 0.7 in units of Mq/Lq , see Fig. 
but within the errors consistent with T derived from the 
longslit data alone. The contours are strongly elongated to- 
wards large T with slightly decreasing M, , as large T in the 
centre can somewhat compensate a lower M,. This degen- 
eracy between M, and T is even stronger for the tiny field 
of view 25mas data alone, where T cannot be constrained 
very well. This emphasizes the importance of including the 
longslit data in the dynamical modelling, which makes the 
degeneracy virtually disappear. 

In order to be able to compare our stellar dynamical Ks 
band T to predictions from population synthesis models, we 
need to scale it to units of M0/L0 by multiplying it with 
the ratio of the band widths: 



Ta- 



y AKs ^ 0.275 yim _ 
AK 0.390 urn ~ 



For a 2 — 3 Gyr old population with a high metallicity 
[ Z/H] = 0.3 - 0.45 (Kuntschner||1998[ |2000[ ) the models 

0.6 (Salpeter 



0. 45 (|Kuntschner||l998 

of |Maraston] ( [19981 12005| ) predicFTA « 0.3 . 
IMF) and Tk ~ 0.2... 0.4 (Kroupa IMF). Thus the mea- 
sured stellar dynamical T is in good agreement with a 
Salpeter IMF but larger than predicted for a Kroupa IMF. 
When taking into account the error in the measured dy- 
namical T and a possible error in the distance assumed for 
Fornax A, then the dynamical T is marginally in agreement 



also with the Kroupa IMF. Cappellari et al. (20061 also find 



that the dynamical T for E/SO galaxies are in agreement or 
larger than the T predicted using a Kroupa IMF. 



5.2 The black- hole mass M. 

5.2.1 M, from the 25mas and the longslit data 

After having constrained the possible T range we now fo- 
cus on the determination of the mass of the central black 
hole. We modelled the 25mas dataset together with the 



longslit data for each quadrant separately and at an in- 
clination of 90°. The resulting values for M, and T with 
the 3(T error bars are listed in Tab. [3] The best-fitting 
black-hole masses and mass-to-light ratios are in the range 
M. = (2.0... 3. 25) X 10* Mo and T = 0.675 ... 0.775 and 



agree very well within 2ct (Axo = X 



Xmin 



6.2) errors. 



5.2.2 M, from the lOOmas and the longslit data 

The same was done using only the lOOmas dataset to- 
gether with the longslit data and the measured PSF. Both 
M, and T are in agreement between the different quad- 
rants, but slightly smaller than for the 25mas data (M. = 
1.0 .. . 2.0 X 10* Mq, T = 0.65 . . . 0.7). Within 2 - 3(j errors 
they still agree with the 25mas results. Note that the sphere 
of influence is well resolved in both cases. 

To figure out whether the uncertainties in the PSF could 
raise a bias in M. , the same models were calculated with 
the 25mas PSF convolved with the kernel from Fig. [2j;. The 
results are virtually identical (see Tab.|3|. Therefore, slight 
uncertainties in the PSF shape do not have a noticeable 
effect on M, or T. 

As a cross-check the 25mas data, binned 4 x 4 to achieve 
the lOOmas spaxel size and convolved with the kernel, were 
also modelled using the convolved 25mas PSF. In this case 
the confidence intervals are much wider and the best-fitting 
values agree within 1 — 2a with both the 25mas and the 
lOOmas results. 

The difference between the two scales could be related 
to the different spatial coverage, although in that case we 
would expect a similar result with larger error bars for the 
25mas models. An increase of T towards the centre could 
also be an explanation. The global T would basically not 
change in that case, and a larger black-hole mass could 
mimic the T gradient. A triaxial structure could also cause 
a systematic difference between the 25mas and the lOOmas 
scale. 

Due to the good overall agreement of the different quad- 
rants this implies that the assumption of axisymmetry is 
justified and therefore the four quadrants of each dataset 
were folded and averaged. The results of the averaged data 
agree with the results of the individual quadrants. 



5.2.3 M, from the longslit and a combination of the 
25mas and the lOOmas data 

In order to take advantage of the high spatial resolution 
of the 25mas dataset and to constrain the orbital distribu- 
tion of the galaxy adequately the 25mas data were com- 
bined with the lOOmas data and this combined dataset was 
modelled together with the longslit data. The spatial region 
covered already by the 25mas data was not considered in 
the lOOmas dataset. The resul ts o f the dynamical modelling 
are shown in Tab. [s] and Fig. 



15 



where Axo is plotted as 
a function of M, and T with error contours for two de- 
grees of freedom. All four quadrants deliver the same results 
within at most 3a errors when all data bins are considered 
(the mean black-hole mass is (M.) = 1.3 x 10* Mq with a 
corresponding rms(M.) — 0.4 x 10* Mq and the mean Ks 
band mass-to-light ratio is (T) = 0.68 with a correspond- 
ing rms(T) = 0.03). When the inner two radial bins are ex- 
cluded the results are almost identical ((A/,) — 1.2 x 10* Mq, 
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Figure 16. Same as Fig. |15| (longslit, lOOmas and 25mas data) 
but for tlie averaged quadrant with (a) or without (b) the bins of 
the central two radii. 
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Figure 17. Fit (solid line) of the best model with M, = 1.5 X 
10* Mq and T = 0.65 to the major and minor axis kinematics 
(points). 



rms(M.) = 0.5 x 10** M©, (T) = 0.70 and rms(T) = 0.02). 
Only the individual Axo contours are wider because of the 
decreased resolution in the centre. 

As the single quadrants agree well with each other it is 
legitimate to just model the LOSVDs folded and averaged 
over all quadrants. This is shown in Fig. [16] for both the 



entire dataset (Fig. 16 1) and with the central bins excluded 
(Fig. 16 a). After the averaging the central velocity peak dis- 
appears, as the LOSVDs of the third and fourth quadrant 
are mirrored. The results therefore are in both cases very 
similar. They also agree very well with the single quadrants. 
The black hole mass of Fornax A, derived from modelling the 
averaged LOSVDs of the combined 25mas and lOOmas SIN- 
FONI dataset and longslit data, is M. = LSlp g'^ x 10* M© 
and the according mass-to-light ratio is T = 0.651q o5^ (3cr 
errors). The fit of the best-fitting model to v, g, hz and /14 
along the major and the minor axis is shown in Fig. |17[ 
Note that the best fit without black hole would be hardly 
distinguishable from the shown fit with black hole (see dis- 
cussion below and Fig. 18 1. The result does not change when 



the inner two radial bins, where the AGN emission distorts 
the CO absorption bands are excluded. Only the error bars 
become somewhat larger (M, = l.S^iJ'j^ x 10* Mq and 
T = 0.725l|]:°25). Note that the la errors derived from a 
smoothed version of the profile agree well with the rms 
of the four quadrant's solution (see above). 

In order to illustrate the significance of the result and 
where the influence of the black hole is largest, Fig. |18| shows 
the differences between the best-fitting model without 



black hole and the best-fitting model with black hole for 
all LOSVDs of the averaged quadrant, analogous to fig. 7 in 
( [Nowak et al.|2007| ). As in the case of NGC 4486a the largest 
black hole signature is found within about 2 spheres of in- 
fluence and in particular along the major axis. For about 
85 per cent of all bins the model with black hole produces 
a fit to the LOSVD better than the model without black 
hole. Adding a black hole improves the fit everywhere and 
not only in the centre, because orbit-based models have a 
lot of freedom and will choose a different orbit distribution 
if no black hole is assumed. Or put in another way: a model 
without a black hole will not just be worse in the very centre 
but will be worse over a relatively large area of the galaxy. 
For the bins with the largest Ax^ along the major axis the 
LOSVD and the fits with and without black hole are shown 
in the left part of Fig. [T8| together with the corresponding 
Axi as a function of line-of-sight velocity. The largest x^ 
differences appear in the high- velocity tails of the LOSVDs. 
The total x^ difference, summed over all LOSVDs, between 
the best model with black hole and the best model with- 
out black hole is Ax^ = 54.7, which corresponds to about 
7.I0". The total x^ values for the models are around 450. 
Together with the number of observables (60 LOSVD bins 
x25 velocity bins) this gives a reduced of ~ 0.3. This is 
a reasonable value as the effective number of observables is 
smaller due to the smoothing ( Gebhardt et al.|2000a K 

For the best-fitting model (T = 0.65) the total stel- 
lar mass within 1 sphere of influence, where the imprint 
of the black hole is strongest, is M. = 1.11 x 10* M©. If 



© 2008 RAS, MNRAS 000,[T|f23| 



The SMBH of Fornax A 15 



0.95 




0.95 




IS 2 25 3 35 

M. (lO^M^;,) 



0.95 
0,9 
0.B5 



0.3 
□.76 - 

0.7 
0.S5 

0,6 
0,55 



0.95 



02 



w % , , . ^ 



0,5 1 1,5 2 2,5 3 3,5 



4,5 




1.5 2 S.5 3 3.5 
M. (10»Mo) 



Figure 15. Models calculated for all four quadrants using a combination of longslit data out to 31 arcsec, lOOmas and 25mas 



SINFONI data. For each quadrant Axq = X 



is plotted as a function of the black-hole mass M, and the Ks band mass- 



to-light ratio T. The coloured regions are the \ . . .ha confidence intervals for two degrees of freedom (Axq = 2.28, 6.20, 12.43, 
19.44 and 28.65). Each model that was calculated is marked as a black dot, the best-fitting model is encircled by a white ring. 
The Axq contours are unsmoothed, which sometimes results in disconnected \a regions due to noise in the models (e.g. in Q3). 



the additional mass of M, — 1.5 x 10 was entirely 

composed of stars, the mass-to-light ratio would increase to 
T = 1.53 (corresponding to Tx = 1.08). This value would 
be typical for an old stellar population (around 8 Gyr for a 
Salpeter IMF, or 13 Gyr for a Kroupa IMF), which has not 
been found in Fornax A ( Goudfrooij et al.|2001 Kuntschner 
2000|. 



Fig. [19] shows the anisotropy profiles for the best-fitting 
model with black hole and the best-fitting model without 
black hole. The models with black hole become tangentially 
anisotropic in the centre (r < 0.2 arcsec), while the models 
without black hole show a certain degree of radial anisotropy 
{(Jr/ot ~ 1.3 where at = [{<Jg + (j^)/2]^^^). This behaviour 
is not surprising, as a central velocity dispersion increase 
can be modelled with either a large black-hole mass or with 
radial orbits. However, the no black hole case is a signifi- 
cantly poorer fit than the best-fitting case with black hole. 
At large radii the tangential bias is decreased with an in- 
creasing black-hole mass, but the anisotropy profile in this 
region is difficult to interpret since we did not include a dark 
halo. 



6 SUMMARY & DISCUSSION 

We have obtained near-IR integral field data with three 
different spatial resolutions 0.085 arcsec for the merger 
remnant Fornax A. Stellar rotation is only detected in the 
outer parts (r > 1.5 arcsec). The stellar velocity dispersion 
slightly decreases from large to smaller radii. In the high- 
est resolution data it shows two peaks in the centre, sep- 
arated by a narrow low-cr region. The average dispersions 
are a = 226 ± 9 km s~^ in an 8 arcsec diameter aperture, 
a = 221 ± 4 km s~^ in a 3.0 arcsec diameter aperture and 
a = 218 ± 8 km s^^ in an 0.8 arcsec diameter aperture. A 
low-luminosity AGN is likely to be present and distort the 
stellar kinematics in the central < 0.06 arcsec, where weak 
[CaVIII] emission in the region of the second CO band- 
head probably induces a velocity peak. The cr-drop can be 
at most partially attributed to this emission line. Either the 
AGN continuum emission or a cold stellar subsystem or both 
could be the cause. 

Near-IR line indices were measured using the high-S/N, 
intermediate resolution data to trace stellar populations. All 
indices show the same trend, a negative gradient and in the 
centre an unresolved depression. The negative gradient of 
the CO index can be interpreted as a metallicity gradient. 
The central drop (~ lA in CO) could be, like the a drop. 
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Table 3. Resulting black-hole masses M, and Ks band mass-to-light ratios T with the according Scr (Axq = 
12.43) errors for all modelled datasets and quadrants. 



Dataset 




Quadrant M, [10^ M0] 


M. _F1 


M.I 


T 






25maM-L^/l 
1 1 1 1 




1 


2.0 


0.3 


3.5 


0.725 


0.625 


0.775 




2 


3.25 


1.5 


4.0 


0.675 


0.625 


0.75 






3 


2.5 


0.7 


3.25 


0.75 


0.675 


0.825 






4 


2.0 


0.3 


4.0 


0.775 


0.675 


0.825 






1 - 4 


1.75 


n 'z 
U. ( 


o.U 


U. ( 


U.Dzo 


Ki. to 






1 - 4 a ( 


^ 3.0 


0.3 


4.5 


0.7 


0.65 


0.75 






]^ 




0.7 


2.5 


0.65 


0.625 


0.7 




2 


1 25 


0.85 


2.75 


0.7 


0.6 


0.7 






Q 


1 


0.3 


2.0 


0.675 


0.65 


0.75 






4 


2.0 


0.5 


3.0 


0.675 


0.625 


0.725 






1 - 4 a 


1.0 


0.3 


1.75 


0.65 


0.6 


0.725 


lOOmas convjTJ- 


-LS 


1 


1.5 


0.5 


2.0 


0.65 


0.65 


0.7 






2 


1 K 


0.7 


2.5 


0.7 


0.65 


0.7 






3 


1.0 


0.3 


1.5 


0.7 


0.65 


0.75 






4 


2.0 


0.5 


2.5 


0.7 


0.65 


0.7 


25mas convj^ -f 


LS 


1 


2.0 


0.0 


4.0 


0.7 


0.65 


0.8 






2 


1.5 


0.3 


4.0 


0.7 


0.65 


0.8 






3 


0.7 


0.3 


3.5 


0.8 


0.7 


0.85 






4 


1.0 


0.3 


2.5 


0.75 


0.7 


0.8 


25mas+100mas 


-fLS 


1 


1.5 


0.3 


2.75 


0.65 


0.625 


0.725 






2 


1.75 


1.0 


3.5 


0.675 


0.625 


0.7 






3 


0.85 


0.3 


2.5 


0.725 


0.65 


0.8 






4 


1.0 


0.5 


1.75 


0.675 


0.625 


0.75 






1-4 a 


1.5 


0.7 


2.25 


0.65 


0.6 


0.725 






1 c 


1.5 


0.3 


3.0 


0.7 


0.625 


0.75 






2 c 


1.75 


0.5 


3.0 


0.675 


0.625 


0.725 






3 c 


0.85 


0.3 


2.5 


0.725 


0.675 


0.8 






4 c 


0.7 


0.3 


2.0 


0.7 


0.65 


0.75 






1 - 4 a c 


1.5 


0.3 


2.75 


0.725 


0.6 


0.75 



Note. When only one SINFONI dataset was included in the modelling, the AI,-T parameter space was 
sampled more coarsely than for the combined SINFONI dataset, with only about half the number of models 
per quadrant. This can sometimes result in the lower or upper limits being identical to the best-fitting values. 



" Lower 3a limit of M. [10^ Mq] 
* Upper 3cr hmit of M. [10** Mq] 

Lower Scr limit of T 

Upper 3(7 limit of T 
^ SINFONI data with a field of field of view of 0.8 X 0.8 arcsec 



f Longslit data from the Siding Spring 2.3m telescope I Saglia et al. 12002 I 
9 The LOSVDs of all four quadrants were averaged and tnen modelled 

Central two radial bins excluded due to AGN contamination 
* SINFONI data with a field of view of 3.0 X 3.0 arcsec together with the measured PSF were modelled 
^ lOOmas SINFONI data were modelled together with the PSF from Fig.[2ji 

^ 25mas SINFONI data, binned to the lOOmas spaxel scale. Data and PSF were then convolved with the 
kernel from 



Fig. [2}; 



caused by the low-luminosity AGN, a subsystem with a dif- 
ferent stellar population, or a combination of the two. 

Using an axisymmetric Schwarzschild code we find a 
black hole with the mass M. = l.Slo g^ 10* ^0 (So" C.L. 
for two degrees of freedom), when the 25mas, the lOOmas 
and the longslit data are included. The same mass with 



slightly larger error bars is found when the innermost bins, 
where the kinematics is distorted by the AGN, are excluded 
from the dynamical modelling. The assumption of axisym- 
metry is justified as consistent results were obtained when 
modelling single quadrants. The mean black-hole mass ob- 
tained from the four single quadrants is {Al,) = 1.3 x 10* Mq 
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(A^-^ = ^2ii ~ Si— 1 (-'(^i noBH ~ -'^iBn) ^.U 25 velocity bins) for all LOSVDs of the averaged quadrant (longslit, lOOmas 
and all 25mas SINFONI data). Bins where the model with black hole fits the LOSVD better are plotted in green the others in 
orange. (b,c) For the radii with the largest positive difference in (a) the LOSVD (open circles with error bars, normalised as 
in |Gebhardt et al.| ( |2000a[ |) and both fits (with black hole, full green line, and without black hole, dashed orange line) are shown 
with the corresponding Ax^ plotted below (d,e). 
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Figure 19. Radial over tangential anisotropy as a function of radius for the best-fitting model with black hole (left) and the 
best-fitting model without black hole (right). The values along the major axis are plotted in black, those along the minor axis in 
cyan and the other colours represent the position angles in between. The vertical dashed line indicates the radius of the sphere of 
influence and the vertical solid line marks the radius out to which the SINFONI data used in the models extend. 



with rms(Af.) = 0.4 x 10* M© and the mean Ks band 
mass-to-light ratio is (T) = 0.68 with rms(T) = 0.03, in 
agreement with the \a limit derived from the analysis of 
the smoothed x^-profile for the averaged data. The PSF of 
the lOOmas data could not be measured very accurately, but 
this does not seem to have noticeable effects. Modelling the 
lOOmas data with the measured or a reconstructed, nois- 
ier PSF (both have about the same FWHM) does not yield 
diflferent results. 



We find a dynamical Ks band mass-to-light ratio of T = 
0.65^0 05^ • When the inner two radial bins are not included it 
is somewhat larger (T = 0.725io'525)- This is in agreement 
with stellar population models assuming a Salpeter IMF, 
but only marginally consistent with a Kroupa IMF. A dark 
halo plays a significant role only outside ~ 31 arcsec, close 
to the effective radius (i?c = 36 arcsec, [Bedregal et al.|2006| ). 
We therefore included neither a dark halo nor longslit data 
at radii larger than 31 arcsec. 
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The black-hole ma ss expected from the M,-a relation 
(Tremaine et al.||2002| is (2.2 ± 0.4) x 10** Mq. Note that 



estimating a^, the luminosity- weighted velocity dispersion 
within Rc, from Fig. [17] gives the same value as a measured 
in the 8 arcsec aperture quoted above. This mass is consis- 
tent with our measurements within the errors. Fornax A is 
only one of two galaxies that underwent a major merger a 
few Gyr ago where a black-hole mass has been measured 
(the other one is CenA, e.g. Neumayer et al. 20071. Both 



seem to be consistent with the M.-cr relation, which holds 
important implications for the growth of the black hole and 
its surrounding bulge. Despite still showing obvious charac- 
teristics of the merging process - like shells and ripples in 
the outer envelope and disordered dust features throughout 
the entire galaxy - the black hole of Fornax A has about the 
mass expected for a normal elliptical. This suggests that 
bulges and black holes grow approximately synchronously 
in the course of a merger, or that the growth is regulated 
during the active phase. 

The black-hole mass of Fornax A is, however, not in 
agreement with the relations between black-hole mass and 
K band luminosity Lk or bulge mass Mbuigc of IMarconi] 



fc Hunt] (|2003|. The total 2MASS Ks band magnitude 
rriKa = 5.587 corresponds to log(Mbuigc/Af0) « 11.4 for 
our Tks = 0.65. For this high bulge mass a black-hole 
mass of ~ 6 X 10* Mq would be expected, a factor ~ 4 
larger than what we measured. Interestingly the situation 
is similar for CenA. With the K band magnitude given in 



Marconi fc Hunt] ( |2003[ ) and Tx=0.72 ( [Silge et al.||2005 l 
the bulge mass is log(Mbuigo/Af0) « 11.0 and a black hole 
with M, 2.2 X 10* M© would thus be expected, a fac- 
tor ~ 5 larger than the mass of 4.5 x 10^ Mq obtained by 



Neumayer et al. (20071. In this context it would be interest- 
ing to reconsider the relationship between M./Afbuigc and 
galaxy age found by Merrifield et al. (20001. Fornax A and 



CenA both have a very small M,/Mbuige and they would be 
among the youngest galaxies in their sample. Although the 
number of merger galaxies with measured black-hole masses 
is fairly small, a correlation between the time of the last 
major merger and M. might start to emerge. 

In order to draw more stringent conclusions on the con- 
nections between merging, bulge growth, black hole growth 
and nuclear activity more dynamical black-hole mass mea- 
surements in merger remnants and in luminous AGN are 
necessary. We were able to show that reliable black hole 
masses can be derived via stellar dynamical modelling even 
if the nuclear source makes dynamical measurements in the 
centre impossible, as long as the spatial resolution is high 
enough that the AGN signature is well within the sphere of 
influence. 
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APPENDIX A: DUST CORRECTION 

The ground-based Ks band SOFI and the H band HST 
NICMOS imaging data were dust-corrected with according 
J band images taken with the same instruments using the 
following method, assuming a standard extinction law: The 
relation between reddening and extinction is 

Ho — He = aEj-H 

for the HST image (analogous for the SOFI image), where 
Ho is the observed, and He is the dust-free image. The ab- 
sorption coefficient a is defined via the extinctions in the 
involved bands: a — Ah,k/{A,j — Ah.k). For the SOFI im- 
age the values Aj = 0.282Av and Ak = 0.112ylv from 



Binney & Merrifield ( 1998 1 were used. The NICMOS image 



was significantly overcorrected with Aj and Ah from |Bin^ 
|ney fc Merrifield] ( |1998[ ) . We found by iteratively changing 
a that the best correction is obtained for a — 0.8. 

The observed images provide the fluxes, so with Eq. 1 
it follows that 

He = Ho - a[{J - H)o - (J - HU 
and 



Ho 



l + a 
Ho 



fa 



Ja 



under the assumption that fj^/fu^ ~ const, (i.e. approx- 
imately no intrinsic stellar population gradient). See also 



Beletsky et al. (20081 
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APPENDIX B: DEPENDENCE OF THE 
RELIABILITY OF THE KINEMATICS ON 
INITIAL AND INTRINSIC PARAMETERS 

In order to find the best value for the smoothing parameter 
a we performed Monte Carlo simulations on a large set of 
model galaxy spectra. These were created from stellar tem- 
plate spectra by convolving them with both Gaussian and 
non-Gaussian LOSVDs and by adding different amounts of 
noise. The results of the simulations might depend, apart 
from the S/N, on the spectral resolution, the LOSVD shape, 
the size of the velocity bins of the LOSVDs 5v, the veloc- 
ity range where the LOSVD is fitted Av, the template it- 
self or other things. Therefore we performed several sets of 
simulations using the best-fitting stellar template for each 
platescale. The spectral range AA — 2.275 . . . 2.349 (xm and 
the number of velocity bins A*' = 29 were identical for all se- 
tups. Sv and Av are given by TV and the velocity dispersion 
cr: Av = 2 X 4.5 X a, Sv = Av/N. We created 100 model 
galaxy spectra for ~ 35 S/N- values between 1 and 170 and 
extracted the kinematics of each spectrum using MPL with 
« 100 different values for a between and 500. We tested 
the dependence on a for different S/N values, platescales, 
templates, a, /13 and /14. 



Bl Dependence on S/N 

The first set of simulations was done using the template star 
HD 12642, observed on the 25mas scale. The stellar spectrum 
was broadened with a LOSVD with the parameters a — 
250 km s~^ and hs — hi = 0. Different amounts of noise 
were added corresponding to a S/N range of 1 . . . 170. The 
kinematic parameters of the resulting spectra were recovered 
using the original spectrum as template star. Fig. |Bl| shows 
the results for three different S/N values. In general, the 
kinematics can be reliably reconstructed for a wide range of 
smoothing parameters [a ~ 0.5 . . . 100 for v and /13, a ~ 
0.5 . . . 50 for a and a « 0.5 . . . 20 for h^). Both a and /14 
change slightly with a. a is biased towards smaller values 
for small a and then starts to increase. /i4 is biased towards 
higher values for small a and continuously decreases and 
is biased towards smaller values for larger a. This trend is 
seen for all S/N. The only difference is in the error bars, 
they increase with decreasing S/N. For small S/N the error 
bars in /14 also increase with decreasing a. As /i4 is more 
sensitive to changes in a than the other parameters, it can 
be used to determine an "optimal" smoothing parameter 
defined as that a where the difference between measured 
and real /14 is smallest. It is shown in Fig. |B2| as a function 
of S/N. For S/N< 25 the curve shows a strong increase and 
then, for S/N< 10, large variations, meaning that it is not 
possible to determine reliable kinematics here. For S/N> 30 
the curve is smooth and decreases very slowly, a sign that 
the kinematics extracted from such data should be reliable. 



B2 Dependence on spectral resolution and 
velocity template 

The spectral resolution of the three platescales is slightly 
different, R ^ 4500 for the 250mas scale, R ^ 5000 for the 
lOOmas scale and R « 6000 for the 25mas scale. This may 
have an effect on the derived kinematics or the optimal a. 



20 40 60 80 100 120 140 160 180 
S/N 

Figure B2. "Optimal" smoothing parameter as a function of 
S/N, derived from the minimum between measured and real 
/i4, the parameter that shows the largest variations with a. In 
order to derive reliable kinematics the S/N of the data should be 
> 30. 



The K5/M0III star HD181109 with a CO equivalent width 
of 12. 9A (see Tab. [l]) has the largest weight in our library 
of velocity templates when fitting to the Fornax A spectra. 
This star was observed in the lOOmas and 250mas scale. 
In the 25mas scale the K5III star IID12642 with a similar 
equivalent width has the largest weight. For each platescale 
spectra with a = 250 km s"'^, /13 — /14 = and different 
amounts of noise as above were created using the according 
stellar template and the kinematic parameters were recov- 
ered with a = 0.1... 500. The results are shown in Fig. 
|B3| for S/N= 140. No significant difference can be detected 
between the three platescales. Likewise the optimal smooth- 
ing parameter is equal for all platescales. Note that this set 
of simulations has been tailored to the Fornax A data and 
cannot be generalised to templates with differing CO equiv- 
alent widths and for velocity dispersions close to or below 
the spectral resolution. 



B3 Dependence on a 

In the next set of simulations the dependence of the ob- 
tained kinematics on the velocity dispersion of the galaxy a 
was tested. The spectrum of HD12642 was broadened with 
LOSVDs with the parameters fta = /14 = and a — 70, 
120, 200 and 250 km s"^ The results are shown in Fig.[B4 



The kinematic parameters for all a values show the same 
trends with a, but with certain differences. With decreasing 
a the LOSVDs become narrower and therefore the scatter 
in v and a decreases. The error bars of /13 and /14 are ap- 
proximately constant with a except for very low dispersions, 
where in addition /14 is biased towards negative values. The 
reason for that could be that the velocity bin width Sv is 
only ~ 20 km s^^ for cr = 70 km s^^ (i.e. smaller than the 
30 km s"'^ px~^ for SINFONI). For low-cr galaxies therefore 
more simulations are necessary to bypass this effect, prefer- 
entially with A*' chosen such that Sv is around 30 km s~^. 



B4 Dependence on and /i4 

With the last set of simulations we tested if the recon- 
structed kinematics depends on /13 or /14. We created model 
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smoothing 



Figure Bl. Kinematic parameters for different S/N (light grey: S/N=40, medium grey: S/N=90, dark grey: S/N=140). The 
kinematic template star HD12642 was used, observed on the 25mas scale. The spectrum was broadened using a LOSVD with 
u = 250 km s~^ and /i3 = /i4 = 0. The reconstructed kinematics is reliable and equal for a wide range of a- values in the range 
of S/Ns of our data (S/N^ 70 for the 250mas and the 25mas and ^ 140 for the lOOmas data on average). As expected only the 
error bars increase with decreasing S/N. 
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galaxy spectra where /13 and hi have the values —0.1 or +0.1 
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B5 



and with a = 250 km s~ . The results are shown in Fig 
in comparison to the case with — = 0. When /is 7^ 
we find deviations from the real values with increasing a 
mainly in v and /13. For simulations with /14 7^ a similar 
behaviour is found for a and hi. 



B5 Results 

Taking all simulations together we conclude that with MPL 
reliable kinematics can be obtained from the CO bandheads 
of the SINFONI data of Fornax A when the chosen a is be- 
tween 1 and 10 and S/N> 30. The S/N of our data is very 
high (on average 70 for the 25mas and the 250mas data and 
140 for the lOOmas data). We are using a = 8 for the 25mas 
and 250mas data and a = 6 for the lOOmas data. These 
results cannot be used unrestrictedly for other datasets, es- 
pecially not for different absorption lines, much lower a close 
to or lower than the spectral resolution, or if the CO equiv- 
alent width is significantly different. 
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Figure B3. Kinematic parameters for the different platescales 
(light grey: 25mas, medium grey: lOOmas, dark grey: 250mas). 
Here the stars HD12642 (25mas) and HD181109 (lOOmas and 
250mas) were used, broadened by LOSVDs with v = km s~^, 
a = 250 km s~^ and = = 0. 
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Figure B5. Kinematic parameters for different and /i4. Here the star HD181109 observed on the lOOmas scale was used. In 
the left column it was broadened by LOSVDs with cr = 250 km s~^ and h'^ = —0.1 (medium grey) or 0.1 (dark grey) and /14 = 0. 
In the right column the stellar spectrum was broadened by LOSVDs with cr = 250 km s~^, /13 = and /14 = —0.1 (medium grey) 
or 0.1 (dark grey). For comparison, the results for /13 = /i4 = are plotted in light grey in both columns. 
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